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Ve  present  a  new  real-time  algorithm  for  the  determination  of  tVie  convex 
hull  of  a  set  of  points  in  the  plane.  Being  real-tim.e,  this  algorithm  has 
the  same  complexity  as  that  of  Preparata  [1],  i.e.  each  new  point  is 
processed  in  time  O(log  N)  v.'here  N  is  the  number  of  points  already  on  the 
convex  hull.  The  method  we  use  is  akin  to  that  of  [1],  but  has  a  simpler 
geometric  structure,  and  has  the  advantage  of  discarding  internal  points 
more  efficiently. 

Both  algorithms  achieve  their  real-time  behaviour  by  means  of  an 
efficient  procedure  to  determine  whether  a  point  P  is  internal  or  external 
to  a  simple  convex  polygon  C,  v/hich  is  the  convex  hvill  of  a  set  S  of 
points  in  the  plane. 

i)  If  P  is  internal  to  C,  then  the  convex  hull  of  S  +  {P}  is  C. 

ii)  Otherwise,  let  vt,  vt'  be  points  on  C  such  that  the  lines  Pvt  and  Pvt 
are  supporting  lines  for  C  (i.e.   C  lies  fully  on  one  side  of  each   line). 
Then  the  convex  hull  of  S+{P}  consists  of  one  of  the  sequences  of  vertices 
of  C  from  vt  to  vt'  (the  sequence  which  is  farther  from  P)   together   with 
the  point  P.  In  what  follows  we  call  vt ,  vt'  the  tangent  points  from  P. 
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to  C.  Otherwise  set  i,  j,  and  k  to  tlie  rnd  points  and   midpoint   of 
the  interval  in  which  the  polarity  changes,  and  go  back  to  3.1 

In  order  to  garantee  that  the  worst-case  behaviour  of  this  algorithm  is 
O(log  N)  we  must  make  sure  that  each  step  does  indeed  divide  each  sequence 
of  vertices  i-'j  and  j-.k  evenly.  This  is  achieved  if  the  hull  C  is 
naintained  as  a  balanced  threaded  tree,  and  if  choose  the  midpoints  1  and 
m  as  the  correct  descondents  of  i,  j  and  k.  The  details  can  be  gleaned 
froni  tVie  following  procedure,  written  in  SETL  [2]. 


PROC  tangent(i , j,k)  ; 

$  The  polarity  of  the  ray  Pv  is  known  to  change  signs  as  v 

$  goes  form  i  to  j  to  k.  The  point  P  is  global  to  this  proc . 

$  The  procedure  returns  the  tangent  point  from  P  v.'hich  is 

$  located  on  the  portion  of  the  hull  from  i  to  k,  if  P  is 

$  external  to  tlic  hull.  If  P  is  internal,  the  procedure  returns 

$  the  SETL  undefined  atom  OM  . 

LOOP  DO 

IF  next(i,j)  and  next(j,k)  THEN 

$  The  three  vertices  i,j,k  are  adjacent. 
RETURN  j  ; 
ENT)  IF  ; 

$  Find  midpoint  between  i  and  j. 

1  :=  IF  Ison(j)  =  i  or  Ison(j)  =  OM  THEN  rson(i)  ?  i 

ELSE  Ison(j)  END  ; 
$  Find  midpoint  between  j  and  k. 

m  :=  IF  rson(j)  =  k  or  rson(j)  =  OM  THEN  Ison(k)  ?  k 

ELSE  rson(j)  END  ; 

pi  :=  polarityd,  P,  1)  ; 
p2  :=  polarityd,  P,  j)  ; 
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p  3  :  =  p  o  1  a  r  i  r  y  (  i ,  V  ,    m )  ; 
pA  :=  polarity (m,  P,  k)  ; 

IF      pi  /=  0  and  pi  /=  p2  THEN  [j,  k]  :-  [1,  j]  ; 

ELSEIv-  p2  /=  p3  THEN  [1,  k]  :=  [1,  m]  ; 

ELSEIF  p4  /=  0  and    p3  /=  pA  THEN  [i,  j]  :-  [j,  m]  ; 

ELSE  return  OM  ;  $  P  is  internal. 

END  IF  ; 

END  LOOP  ; 

END  PROC  ; 

Note  the  following  : 

a)  The  predicate  next(i,j)  can  be  evaluated  in  constant  time  if  the  tree 
that  represents  the  sequence  of  vertices  is  threaded. 

b)  The  special  value  OM  represents  the  undefined  value  (Omega)  and  the 
notation  :  Ison(i)  =  OM  indicates  that  i  has  no  left  descendant  in  the 
tree . 

c)  The  notation  (a  ?  b)  is  equivalent  to  (IF  a  =  OM  THEN  b  ELSE  a) 

d)  the  expression  (polarity  (i,  P,  j)  )  is  evaluated  as 

sign(pix  pj) 
and  requires  two  multiplications  and  a  few  additions. 

e)  The  expression  used  to  determine  the  midpoint  1  can  be  justified  as 
follows : 
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-The  left  descenrfant  of  j  is  the  first  candidate,  ns  it  lies  midway 
between  1  and  j . 

-If  the  left  descendant  of  j  is  i  itself,  then  the  right  descendant  of  i 
lies  nidway  between  them. 

-If  neither  the  left  descendant  of  j  nor  the  ri>;ht  descendant  of  i  exist, 
then  i  and  j  are  adjacent.  Arbitrarily  set  the  midjioint  to  be  i  itself, 
so  that  the  angle  iPl  is  null. 

A  Ronent's  reflection  will   convince  the   reader   that   this   process   is 

correct   regardless   of  the  position  of  i,  j  and  k  in  the  tree,  as  long  as 

the  initial  choice  is  such  that  j  is  midway  betv.'een  i   an   k.   An   initial 

choice  that  has  this  property  is  the  following: 

-  Let  j   be  the  root  of  the  tree. 

-  Let  i   be  the  leftmost  node  of  the  tree,  i.e.  the  first  one  in  order. 

-  Let  k   be   rson( j  ) . 

This  choice  is  also  correct  for  any  circular  permutation  of  i  ,j  ,k. 
Complexity.  * 


The  worst-case  behaviour  of  the  algorithm  is  clearly  0(N  log  N),  given 
that  it  will  take  Log  N  steps  for  the  binary  search  to  determine  that  a 
new  point  P  is  external  to  the  hull  constructed  so  far.  The  average 
behaviour  can  be  seen  to  be  linear  in  N,  for  most  -reasonable  distributions 
of  points  in  the  plane.  This  follows  from  the  results  of  Renyl  [3]  and 
Devroye  [4]  That  state  that  the  expectation  value  of  the  cardinality  C  of 
the  convex  hull  of  N  points  in  the  plane  satisfies: 

E(C)  /  N   ->  0    as  N  ->  <=°. 

so  that  the  complexity  of  our  algorithm  is  governed  by  the  processing  of 
points  that  are  internal  to  the  hull  (This  is  the  case  for  all  convex  hull 
algorithms   based   on   the   throwaway  principle).   The  following  geometric 


argunen*  shows  that   the-   al^^orithra   presonted   above   processes   Internal 
points  in  constant  time  on  the  average. 

After  eacVi  pass  of  our  algorithm,  the  number  of  hull  points  between 
the  end  points  i  and  k  is  halved.  Now  if  the  angle  iPk  is  smaller  than  v, 
then  P  is  immediately  seen  to  be  an  internal  point,  and  tlie  procedure 
terminates.  It  is  only  in  the  case  where  iPk  >  ti  that  a  further  pass  is 
required.  To  estimate  the  probability  that  more  than  n  passes  are 
required,  we  need  to  estimate  the  probability  that  the  vertices  remaining 
under  consideration  subtend  an  angle  of  more  than  it.  This  can  be  done  as 
follows  (Fig. 3)  : 

Let  H  be  the  cardinality  of  the  convex  hull  so  far.  After  n  passes, 
(H  X  2""")  hull  points  remgiin  between  the  end  points  i  and  k.  For  these 
points  to  subtend  an  angle  of  tt  or  greater,  the  remaining  K(  1  -  2  )  hull 
points  must  subtend  and  angle  of  tt  or  less.  If  we  assume  that  the  hull 
points  are  uniformly  distributed  around  P,  then  the  probability  v.'e  seek  is 
equal  to  the  probability  that  given  an  arbitrary  line  through  P,  (H  x  2"") 
points  will  he  in  one  of  the  resulting  half-planes,  and  H(l  -  2"'^)  in  the 
other.  Under  the  assumption  that  the  points  are  uniformly  distributed  in 
the  plane,  the  probability  that  h  of  them  are  on  one  side  is  given  by  the 
binomial  distribution:  (J 

and  the  probabiblity  that  h  or  more  points  are  on  one  side  is  given  by 

x=h  ^   . 

In   successive   passes,   h   takes  the   values   H/2 ,   3H/A,...H(1-   1/2'^). 

Therefore  the  expectation  value  of  the  number  of  passes  needed  is  the   sum 

of   the   probabilities   that   H/2,  3h/4...   points  are  on  one  half-  plane. 
That  is  to  say: 

E(n)  =  1  +  Pn^2  +  P3H/A  +  •  •  • 
m=2     x=h 
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where 


X  =  H(l  -  2~"^) 


To  estimate  this  sun,  we  replace  the  binomial  distribution  by  a  normal 
distribution  with  mean  —  and  variance  -— .  The  sums  above  can  be 
approximated  by  the  integrals; 


log  n 
E(n)  =1+1  j  i/(H/^)  e 

m=2     'Sn 


-Ax^/H 


dx 


Each  term  is  bounded  by: 


A-2/H 


-^x 


./,.2 


p   =  /(h/4tt)  J   e      dx   <    /Ch/Att)  ]   _ 


-^xVH 


m 


m 


dx 


-Ax^/H 


=  /rH76"4ii) 


<    /(H/64tt)  e 


-3H 


l_2-m 
so  that  the  expectation  value  of  the  number  of  passes  is  bounded  by: 

-3H 


3  

E(n)     <  -  +   Log   H    X   /(H/6A7r)e 


<  2 


thus  shov;ing  that  each  internal  point  is  processed  in  constant  time  on  the 
average,  and  therefore  that  our  algorithm  is  linear  in  the  total  number  of 
points . 
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